Critical geometry of two-dimensional passive scalar turbulence 
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Passive scalars advected by a magnetically driven two-dimensional turbulent flow are analyzed 
using methods of statistical topography. The passive tracer concentration is interpreted as the height 
of a random surface and the scaling properties of its contour loops are analyzed. Various exponents 
that describe the loop ensemble are measured and compared to a scaling theory. This leads to a 
geometrical criterion for the intermittency of scalar fluctuations. 
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New experimental techniques which provide complete 
real-space images of fluctuating fields are being devel- 
oped throughout the physical sciences fl]-|j|. They pro- 
vide data on spatial structures inaccessible to methods 
based on system- wide averages, or on measurements that 
probe the system at only a few points. Consequently, 
new measures designed to usefully exploit this wealth of 
information are of great practical importance. Develop- 
ing such measures not only provides the experimentalist 
with new tools with which to analyze and classify corre- 
lations in image data, but may also inspire theories that 
can predict the correlation behavior. 

Scalar turbulence, which deals with dispersion of a 
passive substance (pollutants, temperature, etc.) by a 
turbulent flow, provides an intriguing example of this 
paradigm. When describing the fluctuations of the tracer 
concentration, recently the focus has shifted from the so- 
called structure functions, which characterize the prob- 
ability distribution function for the concentration dif- 
ference between two points, to multi-point correlations 
which better take into account spatial information. This 
has lead to significant advances in turbulence research. 
In particular it has allowed for the identification of the 
source of intermittency of the scalar fluctuations, i.e., the 
breakdown of simple scaling as dictated by dimensional 
analysis [Q. This shift of perspective is nicely exempli- 
fied by a recent study of the Lagrangian evolution of three 
tracer particles ||. It was shown that intermittency in 
the structure functions can be related to the evolution of 
the shape of the triangle defined by the three particles, 
while the overall size of the triangle evolved with a simple 
scaling law. 

Here we develop a geometrical description of scalar 
fields dispersed by two-dimensional turbulence, based on 
their level sets. Unlike earlier studies of the fractal prop- 
erties of the level set , our work focuses on distribu- 
tion functions of level-set loops, in analogy with cluster 
distributions at a critical point. Similar ideas were ex- 
ploited recently by Catrakis and Dimotakis M who intro- 



duced the joint distribution function of shape complexity 
and size to describe level-set islands and lakes obtained 
from images of jet turbulence. We define new measures 
for the random geometry defined by the scalar field and 
propose scaling relations between various exponents that 
characterize these measures. Checking for the validity 
of these scaling relations reveals interesting properties of 
the fluctuations of the scalar field, in particular it pro- 
vides a new, geometrical way of identifying intermittency 
in experimental scalar turbulence data. 

a. Scaling theory of contour loops The scaling 
theory of contour loops of self-affine rough fields (de- 
fined below) was developed in Refs. [|j and ||. Here we 
extend this theory to analyze the scaling properties of 
scalar fields advected by turbulent flows, which are typ- 
ically intermittent, i.e., not self-affine. For this purpose 
we distinguish between scaling relations which are a con- 
sequence of the scale invariance of the loop ensemble only, 
and those which require the self-affine property for the 
concentration field as well. This allows us to test these 
two distinct hypothesis separately against experimental 
data. 

For a random two-dimensional field c(x) , we define the 
contour loop ensemble as a union of level sets for different 
values of c. Each level set consists of contour loops which 
come in many shapes and sizes. The contour loop ensem- 
ble is characterized by the loop correlation function G(r), 
which measures the probability that two points separated 
by r are on the same contour loop, and the joint distri- 
bution function of loop lengths and radii n(s,R). Here 
n(s, R) dsdR is the number of loops, per unit area, which 
pass through a fixed point (say the origin), and whose 
length s e (s, s + ds), and radius R G (R, R + dR). The 
radius of a loop is defined by the side of the smallest 
square that covers it. 

Assuming that the contour loop ensemble is self- 
similar, we can define geometrical exponents which ex- 
press the scaling properties of the two measures, G(r) 
and n(s, R). The loop correlation function, 
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defines the loop correlation exponent Xi, while the joint 
distribution function has a scaling form, 
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and it defines the exponents D and r. Both relations are 
expected to hold for large loops with radii much bigger 
than the lattice spacing (for image data, lattice spacing 
= pixel size). For a self-similar contour loop ensemble 
D is also the fractal dimension of large loops, not too 
be confused with the dimension of the whole level set 
studied previously . 

These three geometrical exponents satisfy a scaling re- 
lation, 



D(3 - t) = 2-2xi 
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which is derived from a sum r ule §. The sum rule 
equates the two ways of calculating the mean loop length 
in a finite region, one using G(r) and the other based on 
n(s,R). This scaling relation is a useful check on the 
validity of the assumption of scale invariance of the loop 
ensemble and in particular Eqs. [l] and |^. 

For a random field which is rough (0 < a < 1) and self 
affinc we can derive a second relation: 



D(t- 1) =2 
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where a is the affine exponent. It is defined by the rela- 
tion c(x) = b~ a c(bx.), where "=" means that the original 
field and its rescaled version are statistically equivalent; 
b > 1 is an arbitrary scale parameter. Eq. (ji|) is a kind 
of hyperscaling relation, and it follows from the fact that 
the number of large loops at a particular scale of obser- 
vation grows with the scale to power a |5). (The usual 
hyperscaling relation, D(t — 1) = 2, say in 2D critical 
percolation, follows from the assumption that at every 
scale the number of large clusters is of order one, i.e., 
a = 0.) 

Furthermore, for self-affine rough fields the three geo- 
metrical exponents, r, D, and xi, depend on the value 
of a only ||. Namely, based on exact results derived in 
the two limiting cases a = and a = 1, x\ = 1/2 inde- 
pendent of a was conjectured J^]; this has been checked 
in numerical simulations |) 10 1. Once this value of xi is 
accepted, formulas for the other two geometrical expo- 
nents, 
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follow from Eq. (g) and Eq. (|). 

We can analyze the level sets of random two- 
dimensional fields based on the above scaling theory. The 
procedure consists of tracing out contour loops of c(x). 



Then for each loop its radius and length are recorded, 
as well as its contribution to the loop correlation func- 
tion. The three geometrical exponents are measured by 
power-law fitting of G(r) (for xi), the average loop length 
as a function of radius (for D), and the number of loops 
whose length is greater than s (for t — 2). The last two 
measures are derived from ro(s, R) by integration: 
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Given that we can reliably extract the three exponents 
(i.e. the data show at least a decade of scaling) we can 
test our assumptions about the loop ensemble by check- 
ing the validity of the various scaling relations. First, 
if the loop ensemble is self similar, Eq. (|J) should hold. 
Furthermore, if it is self-affine, we expect both relations 
in Eq. (||) to hold. We check this by extracting a from 
D and r separately and comparing them. For an inter- 
mittent field we do not expect the two values of a to 
coincide. 

b. Experimental data The contour loop analysis is 
performed on images of a fluorescent dye (passive scalar) 
advected by magnetically forced two-dimensional turbu- 
lent flows (ll]|l^] ■ The flow is generated by running a cur- 
rent through a thin slab of salty water in the presence of 
a perpendicular magnetic field; the Lorentz force on the 
charged ions provides the stirring. The magnetic field is 
generated by an array of permanent magnets placed in 
proximity to the flow. 

With a suitable arrangement of magnets and choice of 
driving current a turbulent velocity field is set up either 
in the regime of a direct enstrophy cascade |ll|| or an 
inverse energy cascade 0|. In both cases the advection 
of fluorescent dye was monitored using a CCD camera. 
The data we analyze below consists of two sets, one for 
each of the two flow regimes, of 10 images of the tracer 
field. These images were selected from the time interval 
during which stationarity of the tracer fluctuations was 
achieved. In this time interval the energy dissipation was 
observed to be roughly time-independent |nj . 

c. Results of the loop analysis The contour loop 
analysis for the two types of forcing was performed by 
measuring the loop correlation function, the average loop 
length as a function of the loop radius, and the distribu- 
tion of loop lengths from the images provided by P. Tabel- 
ing's group. Each image consists of a 512 x 512 array of 
integers which represent the gray scale values of the in- 
tensity of the fluorescent dye, which is proportional to its 
concentration c(x) . The unit of length we use throughout 
is the lattice spacing which physically corresponds to the 
size of one pixel. For each image 10 4 points on the dual 
lattice were chosen at random and through each point 
a contour loop was constructed following the algorithm 
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described in ||. For each loop we measured the length 
and radius, and its contribution to the loop correlation 
function. 

To check the stationarity condition for the data from 
the two flow regimes we performed the loop analysis on 



each image separately. We found no discernible variation 
in time for the three measures G(r), (s) (R), and N > (s). 
Therefore we combined the loop data from all 10 images 
to improve our statistics. 
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FIG. 1 The top row gives the level set and loop analysis of the scalar data in the direct cascade while the bottom row 
figures corresponds to scalar data obtained from the inverse cascade. All lengths (r,R, and s) are reported in units of 
one pixel; G(r) is not normalized. 



i. Direct cascade - The loop correlation function shows 
the 1/r (i.e. xi — 1/2) behavior characteristic of rough 
fields for separations r < 100 (see Fig. 1); in terms of 
the geometrical exponent xi we find from a linear least 
squares fit a value of xi — 0.495(10). The rather large 
(compared to the inverse cascade data) crossover region 
for 100 < r < 300 which precedes the fall off at r > 400 
(due to system size) might be indicative of an intermedi- 
ate length scale at which the scaling of the loop ensemble 
changes. This is confirmed by the cumulative loop length 
distribution, N > (s). 

The graph of the mean loop length as a function of 
the radius, Fig. 1, shows a scaling range for loop radii 
3 < R < 300. It leads to a fractal dimension for contour 
loops of D = 1.20(3). The value of the exponent is esti- 
mated by comparing results from linear least-square fits 
to data within several different subintervals of the scaling 
range. 

The cumulative loop length distribution graph, Fig. 1, 



shows scaling for 10 < s < 200 after which there is a 
crossover to what looks like a power law with a smaller 
(more negative) exponent. The fall off at lengths s > 
7000 occurs as the system size intervenes. The crossover 
occurs at length scales comparable to the effective dissi- 
pation scale which was identified from the scalar power 
spectrum in Ref. p^ . The exponent r — 2 = 0.19(2) is ex- 
tracted from the scaling range, using the same procedure 
(sub- interval fits) as for Df. 

Using the scaling relations Eq. (JsJ) , we extract the value 
of the affine exponent a = 0.60(6) from the measured D, 
and a = 0.53(8) from t — 2. The fact that they match 
within error bars indicates that the loop data is consistent 
with simple scaling of the tracer field below the effective 
dissipation scale. 

Above the dissipation scale Batchelor (logarithmic) 
scaling is expected to hold |ll],|l3|]. This implies a = 
and is not well supported by the loop analysis. Namely, 
we would expect no crossover for xi which is universal, 
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and (from Eq. (J5|) ) a crossover of D to higher and r — 2 
to a lower value. This is qualitatively consistent with the 
loop data (accept for x{) but the actual numbers are off. 
We conclude that if the Batchelor regime is reached it 
occurs for a narrow range of scales which is not resolved 
by the loop analysis. 

Additional checks on the critical nature of the loop en- 
semble is provided by measuring the distribution of loop 
radii, N > (R). N > (R) is the number of loops whose ra- 
dius is larger than R and N>{R) ~ R- D( - T -^ follows 
from Eq. (||) . In the AT> (R) plot we again see a crossover 
occurring at R « 60, which from the (s) (R) plot in Fig. 1 
corresponds to s ~ 300, consistent with the N > (s) plot. 
For loops with 3 < R < 40 we measure D(r-2) = 0.22(1) 
by fitting subintervals to a power law, while from the pre- 
vious measurements of D and r — 2, D(t — 2) = 0.23(3) 
follows. 

ii. Inverse cascade - Loop analysis of the inverse energy 
cascade images reveals a very different picture. The scal- 
ing is much better as no crossover scale is detected. The 
geometrical exponents can be extracted reliably from the 
data and they do not satisfy the scaling relations Eq. (||) . 

The loop correlation function decays with distance as a 
power law, for r < 100 (see Fig. 1); the measured value of 
the exponent xi = 0.485(5) is consistent with the univer- 
sal value xi = 1/2. The (s) (R) plot shows a scaling range 
for 10 < R < 100, from which we extract D = 1.30(1) 
(see Fig. 1). The N > (s) plot, Fig. 1, has a rather re- 
markable scaling range for loop lengths 10 < s < 3000 
from which the geometrical exponent r — 2 = 0.273(4) 
is obtained. Assuming the validity of the affine scaling 
relations Eq. (pj) the exponent a = 0.40(2) follows from 
the value of D, while a = 0.25(2) is the value implied by 
the measured t — 2. Clearly, these two values are differ- 
ent, thus signaling the breakdown of simple scaling. We 
conclude that the scalar field is intermittent. 

To check that the loop ensemble is described by a joint 
distribution function that has a scaling form, Eq. (2), 
we check for the validity of the scaling relation Eq. (2). 
From the measured values of D and r — 2 we have 
D(t — 3) = 0.95(1) while we expect 2 — 2x\ = 1 from 
the universal value xi — 1/2; 2 — 2xi = 1.03(1) is the 
value that follows from the G{r) measure. Furthermore 
we measure the combination D(t — 2) from N > (R). We 
find D(t - 2) = 0.344(4) (scaling range 3 < R < 100), 
while the measures used above yield D(t — 2) = 0.355(8). 
Taken together these results give strong support to the 
assumed critical nature of the contour loop ensemble of 
the scalar data in the inverse energy cascade regime. 

In conclusion we have shown that a scaling analysis of 
contour loops of a passive scalar field advected by turbu- 
lent flow is a useful geometrical measure for studying its 



fluctuations. The loop analysis successfully distinguishes 
between different flow regimes and indicates the break- 
down of simple scaling in the case of the inverse cascade. 
Even though the scalar field is intermittent in this case, 
the loop ensemble is critical and characterized by two in- 
dependent geometrical exponents (D and r). The direct 
cascade data reveals an intermediate length scale which 
can be identified with an effective dissipation length. Be- 
low this dissipation scale simple scaling of the scalar field 
is favored by the loop analysis, while the data above this 
scale is not indicative of a Batchelor regime with loga- 
rithmic structure functions. 

The loop analysis described here has an advantage over 
structure-function measurements as the geometrical ex- 
ponents are sensitive to higher moments of the concen- 
tration difference and they capture the spatial structure 
of the fluctuations more fully by incorporating connect- 
edness properties. As three dimensional data becomes 
more available extensions of these ideas to contour sur- 
faces will be a promising direction for future research. 
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